LATE-TIME BEHAVIOUR OF THE TILTED BIANCHI TYPE VI 1/9 MODELS 
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Abstract. We study tilted perfect fluid cosmological models with a constant equation of state parameter 
in spatially homogeneous models of Bianchi type VI_i/g using dynamical systems methods and numerical 
simulations. We study models with and without vorticity, with an emphasis on their future asymptotic 
evolution. We show that for models with vorticity there exists, in a small region of parameter space, a 
closed curve acting as the attractor. 



1. Introduction 

In recent papers the tilted Bianchi models [H [21 El HI 13 El 13 IH1 El LlOl L^] have been studied using dynamical 
systems methods (T2J [T3] . In this paper we will study the general tilteqj perfect fluid Bianchi model of type 
VI_i/9, with a constant linear equation of state parameter 7 satisfying the causality conditions (i.e., no 
super luminal speed of sound), which was not studied in detail in [S]. The tilted models have shown a wide 
variety of phenomena including, for example, the existence of closed curves [8] and incompleteness of the 
fluid congruence [15l I16j . 

The general irrotational perfect fluid type VT_i/g model is the most general of the irrotational Bianchi 
models. In fact, the irrotational type VI_i/g model has a 6-dimensional state space; just one dimension less 
than that of of the most general tilted type VI_i/g model. We will show that, regardless of whether we 
consider the general type VI_!/ 9 model with or without vorticity, the important late time asymptotes are 
the non-tilted Collins type VI_!/ 9 solution (for 2/3 < 7 < 10/9), 

9 9 9 9 2 ( 4 ~ 3 ~0 2sa 9 _4_ 4sx 9 

ds = -dt + t dx + t 3 ~> e 3 ~< dy + t 3 -ie 3 ~< dz , 

where s — y(2 — 7) (37 — 2), or the Collinson- French (or Robinson- Trautman) vacuum solution (for 10/9 < 
7 < 2), 

ds 2 = -dt 2 + t 2 dx 2 + ti (V^dy + ^tidx^j + tie^dz 2 . 

The tilt of the fluid can be asymptotically non-tilted, intermediate or extreme depending on the value of 7. 
The Collinson-French solution is a Petrov type III vacuum solution, and hence it has the peculiar feature 
that all the curvature invariants of zeroth order vanistQ. 

The spatial hypersurfaces in a Bianchi type VI_i/g cosmology are defined as the orbit of the action of the 
Bianchi type VL_i/g Lie group. This Lie group can be equipped with the left-invariant one-forms: 

= dx, Z> 2 = e~ bx dy, £ 3 = e 2bx dz, 

where b is a constant. The exceptional feature of this model does not relate to the group itself (for which 
there is nothing extraordinary) but arises from the Einstein field equations. For this particular model one 
of the Einstein constraint equations vanishes exactly and hence the vacuum case allows for an extra shear 
degree of freedom. This extra shear is present in the Collinson-French vacuum solution and therefore this 
metric has no analogue in the other Bianchi vacuum cosmologies. 

The aim of this paper is to fill one of the gaps in the analysis of the general tilted type VI^ models 
presented in [9]. The h = —1/9 case is special (as explained above) and requires a separate analysis, which 
has only been done in part previously. In particular, we will show that closed periodic orbits exist for the 
general vortic model. However, for the irrotational models, no such closed curve exists. Moreover, we will 



Date: February 1, 2008. 

^^The non-tilted Bianchi model of type VI_]/g was studied in 1 14| . 

2 However, considering differential invariants they will not vanish; for example, C a p 1 s-^C a P~ ,s ''^ ^ 0. 
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also complete a centre manifold analysis appearing in one 6-dimensional invariant subspace describing vortic 
Bianchi type VI_i/g models. All other attractors will be included for completeness. 



2. Equations of motion 
2.1. The orthonormal frame approach. The line-element of a Bianchi cosmology can be written 

(2.1) ds 2 = -dt 2 + 6 ab LO a u b , 

where t is the co-moving cosmological time. The one-forms ui a are left-invariant one-forms on the hyper- 
surfaces spanned by the group orbits and can be related to the left-invariant one-forms, a; 1 , given above by 
u a = e a t (t)u. 

The geometric (or normal) congruence, n M , is given by n = d/dt. It is also useful to define and the Hubble 
scalar and the shear associated with the congruence n^: 



(2.2) 



H = -n* 
3 



' [If 



The matter variables are chosen to be the energy density, (i, and the tilt-velocity, v a , which is defined as the 
3- velocity of the fluid with respect to the geometric (or normal) congruence, nf*. The equations of motion 
can now be written down in terms of the Hubble scalar, H ; the shear, a a b] the curvature variables n ab and 
a c ; and the matter variables /i and v a . 

Expansion-normalised variables are introduced (by dividing the variables with the appropriate powers of 
H) . The papers [7J QT] contain all the details regarding the determination of the evolution equations for 
the tilted cosmological models under consideration. Here, we shall adopt the so-called iV-gauge in which 
the function N x is purely imaginary; this is realised by the choice <fi' — \/3AE_, where A is defined by 
N = Alm(Nx). The evolution equation for TV can then be replaced by an evolution equation for A, which 
ensures a closed system of equations. We will also adopt the dimensionless time parameter r, which is related 
to the cosmological time t via dt/dr = (1/H), where H is the Hubble scalar. 

Using expansion-normalised variables, the equations of motion are (see [Tllllj for the complete derivation 
of the equations): 



(2.3) S' + = ( g -2)S + + 3(E 2 2 + E 2 3)-27V 2 + ^-(~2 W 2 + « 2 2 + t- 3 2 ) 

(2.4) E'_ = (q-2-2v^S 23 A)S_ + v^(I] 2 2 -S 2 3) + 2A7V+^^(v 2 -w 2 ) 

(2.5) S' 12 = (g-2-3S+- V3£-)£i2- V3 (S 23 + S_A) S 13 + ^J^ Wl ?; 2 

(2.6) E' 13 = ( 9 -2-3S + + V3S_)s 13 -v / 3(S 23 ~E_A)I]i 2 + ^|^wi«3 

(2.7) S 23 = ( (Z -2)S 2 3-2v / 3A^ 2 A + 2V3AE 2 _+2V3Si 2 Ei3 + ^^W2W3 

G+ 

(2.8) N' = (q + 2S+ + 2V3S 23 a) N 

(2.9) A' = 2\/3E 2 3 (1 - A 2 ) 

(2.10) A' = (q + 2E + )A. 
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The equations for the fluid are: 

= ^{2<Z - (3 7 ~ 2) + 2 1 Av 1 + [2g( 7 - 1) - (2- 7 ) - 7 <S] V 2 } 

= (T + 2£+) Wl - 2V3S13W3 - 2v / 3Si2« 2 - A (»| + w|) - V3JV (w 
= (T - £ + - v/3S- 



(2.11) 


n' 


(2.12) 


v[ 


(2.13) 


v' 2 


(2.14) 


v' 3 


(2.15) 


V 


where 





v 2 -Vz (S 23 + £-A) u 3 + \/3AiVui«3 + (A + 
T - £+ + V3£-) u 3 - v 7 ^ (£ 23 - £- A) v 2 - v^ViVviVa + (a - 

y(i-i/ 2 ) 



ViV 2 

Vzn) v lV3 



i - ( 7 - i)y 2 



[(37 - 4) - 2( 7 - l)A Wl - 5] , 



q = 

£ 2 = 

s = 

V 2 = 

T = 

G+ = 



2£ 2 

^2 



l(3 7 -2) + (2- 7 )F 2 



2 

£ 2 



n 



e 2 

^13 



^23 



i + ( 7 -i)y 2 

£^ + £^ + £ 2 2 

£af>C C , C C a = 1, 
2 i 2 i 2 

[(3 7 - 4) - 2( 7 - (1 - F 2 ) + (2 - 7 )F 2 S 

1 - ( 7 - 1)V 2 

l + ( 7 -l)^ 2 . 



These variables are subject to the constraints 



(2.16) 
(2-17) 

(2-18) 

(2-19) 
(2.20) 



= £ 



■ A 2 + A^ 2 + 



= 2E+A + 2E_7V- 



7 rjfi 



= - Si 2 (JV + V3A) + Ei 3 AJV 

= [Ei 3 (JV- V3A) + E 12 AiV 
= 3A 2 - (1 - A 2 ) iV 2 . 



7 ^w 2 

7^"3 



The parameter 7 will be assumed to satisfy 7 € (0, 2). The generalized Friedmann equation (|2. 16|) yields an 
expression which effectively defines the energy density Q. We will assume that this energy density is non- 
negative: f2 > 0. Therefore, the state vector can be considered X = [£ + , £_, £i 2 , £13, £ 23 , N, A, A, v\, v 2 , v 3 ] 
modulo the constraint equations (|2. lT[) - (|2.20p . Thus the dimension of the physical state space is seven (and 
hence of equal generality to the other models of type VI^ for a given value of h). Additional details are 
presented in [7]. 

The dynamical system is invariant under the following discrete symmetries : 



[£+, £_, £12, £13, £ 23 , N, A, A, v 1 ,v 2 ,v 3 ] h-> [£+, £_, £i 2) £13, £23, -N, A, -A, -v 1} -v 2 , -v 3 ] 
[£+, £_, £12, £13, £23, N, A, A, vi,v 2 ,v 3 
[£+,£_, £12, £13, £23, N, A, A, Vi,v 2 ,v 3 
[£+, £_, £12, £13, £23, N, A, A, Vi,v 2 ,v 3 



[£+, -£_, £13, £12, £23, -N, A, A, v%,v 3 , v 2 ] 
[£+,£_, ±£12, T£i3, -£23, N, -A, A, Vi,±v 2 , Tv 3 ] 
[£+, £-, -£12, -£13, £23, N, A, A, vi,-v 2 , -v 3 ] 



These discrete symmetries imply that without loss of generality we can restrict the variables A > and 
N >0, since the dynamics in the other regions can be obtained by simply applying one or more of the maps 
above. The third and fourth symmetries listed imply that one can add additional constraints on the variables 
£12, £i3,i>2 or d 3 . For the type VL1/9 we can use </>4 to restrict w 2 to be non-negative (since v 2 = implies 
w 3 = 0, see below); hence, we will assume w 2 > 0. There is no natural way to restrict any of the remaining 
variables using the symmetry ^3 , and so we will not do so here. 
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2.2. Invariant sets. For the case h = — 1/9 we define (N a bV a v h is identically zero for h = —1/9) 
(2.21) 3= [A(E2 2 + S2 3 ) + 2S 12 E 13 ] . 

In this analysis we will be concerned with the following invariant sets: 

(1) T(VI_i/ 9 ): The general tilted type VI_ 1/9 model: |A| < 1. 

(2) Ti(VI_ 1/9 ): A one-tilted type VI_ 1/9 model: |A| < 1, v 2 = v 3 = S i2 = S i3 = 0. 

(3) Ti. (VI_i/g): A one-tilted diagonal type VI_i/g model: v 2 = v 3 — £i 2 = £13 = £23 = A = 0. 

(4) N ± {VI- 1/9 ): A class of tilted type VL_i/g models: |A| < 1, D = 0. 

(5) T 2 (V I-i/g): A two-tilted type VI_i/ 9 model. This is the fixed-point-set of (j) 3 and is given by 
v 3 = E13 = £ 23 = A = 0. 

(6) T 2 ~ (V I_i/ 9 ): A one-tilted irrotational type VI_!/ 9 model. This is the fixed-point-set of 4> 3 and is 
given by v 2 — v 3 = S i2 = £ 23 = A = 0. 

(7) B(VI_ 1/9 ): Non-tilted type VI_ 1/9 : |A| < 1, Vl = v 2 = v 3 = £ 12 = S 13 = 0. 

(8) B (VI_i/ 9 ): A class of diagonal non-tilted type VI_ 1 / 9 models (n a a =0): V = Si 2 = S13 = S 23 = 
A = 

(9) T(II): The general type II model: A = ±1, A = 0. 

(10) B{I): Type I: N = A = V = 0. 

(11) dT(I): "Tilted" vacuum type I: fi = N = A = 0. 

Regarding A r± (V r /_ 1 / 9 ), to verify that this is indeed an invariant subspace, we calculate D 1 : 



D' = 2(q-2- 3S+ - V3AE 23 + ZAvi^J D; 



hence, D — defines an invariant subspace. Note also that this invariant set is not a manifold; it is similar 
to the light-cone in 2-dimensional Minkowski space. Therefore, N ± (VI- 1/9 ) - Ti(VI- 1/9 ) consists of 4 
disconnected pieces. By the symmetry cj>^, these are actually only two inequivalent pieces. Here, we choose 
N ± (VI_ 1/9 ) such that 

T+(VI_ 1/9 ) c N+(VI_ 1/9 ), T 2 -(VI_ 1/9 ) c N-(VI_ 1/9 ) 

Since N+(VI_ 1/9 ) n N-(VI_ 1/9 ) = T 1 (VI_ 1/9 ) 1 both N+(VI^ 1/9 ) and N-(VI_ 1/9 ) arc invariant sets. 
We note that the closure of the set T(VI_i/ 9 ) is given by 



(2.22) T(VI_ 1/9 ) = T(VI_ 1/9 )UT(II)uB(I)UdT(I). 

Since the boundaries may play an important role in the evolution of the dynamical system we must consider 
all of the sets in the decomposition (|2.22[) . 

Let us consider the constraint equations (|2.18p and (|2.19[) as a linear map 

L: (S 12 ,Si3) h-> (V2,v 3 )/G+, 

where L is considered as given in terms of A, N, A and f2. For h ^ —1/9, det(L) ^ and the image of L is 
2-dimensional. However, for h = —1/9, det(L) = and the image of L is 1-dimcnsional; hence, in this sense 
the constraint equations are degenerate. This implies that (w 2 ,f 3 ) has to be restricted to a 1-dimensional 
submanifold. We will therefore say that the general type VI_i/ 9 model only allows for 2 tilt degrees of 
freedom. In particular, we can solve for v 3 and obtain 

_ \v 2 

It is illustrative to consider the eigenvectors of the map L: 

xo = (- TT -A = =, l) , Lx = 0, 

1 ' a! „ ) > Lx « = ax a, a > 0. 

l + Vl-A 2 / 

For each of these eigenvectors, D = and hence, alternatively, we can define iV ± (y/_ 1 / 9 ) when (Ei 2 ,Ei 3 ) 
is proportional to one of these eigenvectors. More specifically, for N + (VI-i/ 9 ), (Ei 2 ,Si 3 ) cx x a , while for 
N-(yi- 1/9 ), (S 12 ,E 13 )cxx . 
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Table 1. The dimensions of the invariant sets for the Bianchi type VI-ug model. The 
right-most column indicates the specialization in terms of the Gi cosmologies. Here, Diag 
means diagonal, HO means hypersurface orthogonal, and OT means orthogonally transitive. 
The stars indicate the exceptional cases where the models aquire an addition shear degree 
of freedom. 



Dim 


Invariant set 


No of Tilt 


Gi action 


2 


B (VI_ 1/9 ) 





Diag 


3 


Ti, (W_i /9 ) 




1 


Diag 




Bo(vr_ 1/g ) 


■k 





HO 


4 


B(VI_ 1/9 ) 







OT 




T?(VI^ 1/9 ) 




2 


HO 




T 2 -(VI_ 1/9 ) 


* 


1 


HO 


5 


Ti(VI_ 1/9 ) 




1 


OT 




B{vr_ 1/9 ) 


-*- 





Gen G 2 


6 


N+{VI_ 1/9 ) 




2 


Gen G 2 




N-(VI_ 1/9 ) 


* 


1 


Gen G 2 


7 


T(VI_ 1/9 ) 


2 


Gen G 2 



Note also that, due to the presence of the eigenvector with zero eigenvalue x , (^2,^3) = does not 
necessarily mean that (£12, £13) is zero. In particular, for the non-tilted models this implies that we may 
have an additional shear degree of freedom; these models have usually been called the exceptional case and 
are denoted by an asterisk; e.g., B 9 (VI*_ l ^) and B(VI*_ 1 ^ g ). We also note that for the tilted models, there 
is an exceptional case of the one-tilted models Ti(VI_i/ 9 ) which could be denoted T^VT^ , g ). However, 
Tx(V r J* 1 , 9 ) = N~(VI-i/ 9 ) as explained above. Therefore, we keep the notation N~(VI-i/ 9 ). Similarly, we 
have T lt0 (VI*_ 1/9 ) = T a -(VJ_ 1/9 ). 

2.3. Fluid Vorticity. The various invariant subspaces can also be categorised in terms of the (-f/fluid- 
normalised, where iffluid = fluid vorticity, W a . The vorticity of the fluid for the type VI_i/ 9 

models is given by: 

(2.23) W a = ^-(N ab + s abc A c )v\ W = 0, 

where 

1 - \{V 2 + V 2 S + 2A a v a ) 



B = 



[1 - {j-i)v 2 }Vi-v 2 

For the invariant sets: 

(1) T{VI_ 1/9 ): W° = W 1 = 0, most general vortic type VI_ 1/9 . 

(2) N+{VI-i/ 9 ): W° ^W 1 = 0. 

(3) N-(VI- 1/9 ): W° = W a = 0, non-vortic. 

(4) T+(y/_ 1/9 ): W° = W 1 = W 2 = 0. 

(5) r 2 ~(V/_i/ 9 ): W° — W a = 0, non-vortic. 

(6) Ti(F7_i/ 9 ): W° — W a = 0, non-vortic. 

(7) B(VI_ 1/9 ): W° = W a = 0, non-tilted and non-vortic. 

We note that the most general non-vortic model is of dimension 6. Hence, since the non-vortic type VI^ 
and Vllft, models - regarding h as fixed - are of dimension 5, the type VT_i/g model is the most general 
non-vortic model of all Bianchi models. 

In general we can also solve for the vorticity component W 2 : 



W 2 = W s 

This follows from the constraint equations and eq. (|2.23D . 



A 
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3. Qualitative behaviour 

3.1. Monotone functions. There are a number of monotone functions in the state space of interest. For 
< 7 < 6/7, there exists a monotonically increasing function Z\ defined by 

(l-7 2 )3(2-7) 

(3-1) Zi = art~\ a = K Gl _\ vi , 

Z[ = [ 2 (l- 7 )g + (2-7)+ 7 5]Z 1 . 
This function can be used to show [7]: 

Theorem 3.1. For < 7 < 6/7, aZ/ tafted Bianchi models (with Q > 0, V < 1) of solvable type are 
asymptotically non-tilted at late times. 

Corollary 3.2 (Cosmic no-hair). For Q > 0, V < 1, and < 7 < 2/3 we have that 

lim n = 1, lim 7 = 0. 



Moreover, the following function is a monotone function in T(V/_ 1 yg): 

A A N 2 GlD 2 



(3.2) Z 2 - 

V 1 (l-T/ 2 )5(2-7)f}5' 

(3.3) = (57-6)(3-2A«!)Z 2 , 

where is defined in cq. (|2.21[) . This function is monotonically decreasing for 7 < 6/5 and monotonically 
increasing for 6/5 < 7. 

We note that in the subspace N~ (V T_i/g) we have the monotone function: 

(3.4) Z 3 - 



A 2 G+(1 -1/2)^(2-7)' 

(3.5) = -(2- 7 )(3-2Au 1 )Z 3 . 

This function is monotonically decreasing in N~ (V I-i/q). 
The monotonic function Z3 immediately implies: 

Theorem 3.3 (Future behaviour in N-(VI- 1/Q )). For 2/3 < j < 2, fl > 0, A > 0, v 2 < I, v 2 = v 3 = 
have that: 

either lim f2 = 0, or lim 7 = 0. 



r — >oo 



This implies that all irrotational type VI_!/g universes are either asymptotically vacuum or non-tilted at 
late times. 



3.2. Equilibrium points. 

3.2.1. B(I): Equilibrium points of Bianchi type I. 

(1) 1(1): S + = E_ = E12 = £13 = £23 = A = N= V = and = 1. Here, |A| < 1 and is an 
unphysical parameter. This represents the flat Friedman-Lemaitre model. 

The remaining equilibrium points are all in dT(I). 

3.2.2. T(II): Equilibrium points of Bianchi type II. All of the tilted equilibrium points come in pairs. These 
represent identical solutions (they differ by a frame rotation); however, since their embeddings in the full 
state space are inequivalcnt, two of their eigenvalues are different. All of these equilibrium points have an 
unstable direction with eigenvalue — 2v / 3£23 corresponding to the variable A. These equilibrium points are 
given in [5]. 
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3.2.3. T(VI-x/g): Equilibrium points of Bianchi type VI-i/g. 

(1) C(-l/9): Collins perfect fluid solution, 2/3 < 7 < 5/3 

Eia = E13 = £ 23 = A = V = 0, £+ = -f(37 ~ 2), £_ = ^f(3 7 - 2), N 2 = ^(3 7 - 2) (2 - 7), 
A = N/V3, ft = |(5 - 37). This equilibrium point is in B(V7_i/ 9 ). 

(2) 7£+(-l/9): The Apostolopoulos h = -1/9 solution [17], 4/3 < 7 < 3/2 

This is a vortic solution lying in the invariant subspace T + (V The solution is given in terms 
of the expansion-normalised variables in |llj with h = —1/9, k = 1/3. 

(3) Bianchi type VI_i/g vacuum plane waves. All of these solutions have 

ft = £12 = E13 = £23 = 0, £_ = N = v / -E+(l + £+), A = (1 + £+), -!<£+< -3/4, |A| < 1. 



It is avantageous to introduce r = \fl — A 2 , which implies that we can write 

z+ = -^rr^' o<r<i. 

3 + r z 

We will also define p by 

p = V | + ?;|. 

The equilibrium points are then determined by the tilt velocities: 

(a) £(—1/9): v\ = v 2 = W3 = 0. These represent 'non-tilted' plane waves and lie in B(VI_x/g). 

(b) £(-1/9): V! = 37(3+ 2 ;' ) ~^ ( 1 9 ) +2r2) , v 2 = v 3 = 0, ^±rji < 7 < 2 . These represent 'intermediately 
tilted' plane waves and lie in Ti(V7_x/g). 

(c) £±(—1/9): v% — ±1, «2 = V3 = 0. These represent 'extremely tilted' plane waves and lie in 
Ti(W_ 1/9 ). 

(d) F+{h): Here, (9 + 7r 2 )/[3(3 + r 2 )] < 7 < 3(3 + r 2 )/(9 - r 2 ) and 



37(3 + r 2 ) - (9 + 7r 2 ) 
2r 2 (3- 7 ) : 
(9 + r 2 ) [5 - 3 7 ] [3 7 (3 + r 2 ) - (9 + 7r 2 )] 



Vl = 2^3^ ' V *- V * =pr 



8r 4 (3 - 7 ) 2 

These represent 'intermediately tilted' plane waves and lie in N + (VI_ 1/g ) (for A = they lie 
in T+(VI_ 1/9 )). 
(e) £+(-1/9), < 7 < 2: 

4r 2 

vi = ~w^r vl ~ vl = pr 

2 (9 + r 2 )(9-7r 2 ) 
P = !-«!= 9(3 - ,2)2 

These represent 'extremely tilted' plane waves and lie in AT + (W_i/g). For A = these equilib- 
rium points lie in T+(V7_i/9), and due to the special importance for the late-time behaviour, 

we will denote £+(-1/9) = £+(-1/9) . Therefore, 

A— 

2 



£+(-1/9): «! = --, «a--g-. 
(4) C7 7 : The Collinson-French (Robinson- Trautmann) solution is given by: 



£12 = £23 = ft = A = 0. 

There are the following equilibrium points associated with the Collinson-French solution: 
(a) CT a : Vi = v 2 = 0, < 7 < 2. 



/v\ rir %/6(3 7 -4) v / 5(3 7 -4)(3-2 7 ) 4 

(c) C^ 2 : «i = 4^,«2 = 0,M^ <7< M±j/I. 
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(d) £CT±: Vx = ±1, v 2 = 0, < 7 < 2. 

These equilibrium points, and their stability, were studied in [TJ. 
(5) W: Wainwright 7 = 10/9 solution: 



+ 3' 



1 



3V3' 

S12 = S 23 = A = V = 0. 



< £13 < 



15 



6 



54E? 3 , 



n = 5 - 
9 



3£ 



13; 



4. Late-time behaviour 
The late-time behaviour of models with < 7 < 2/3 is determined by (the Cosmic no-hair) Corollary 3.2. 

4.1. The invariant subspace N + (VI-i/g). Here, we have the following late-time attractors: 

• 2/3 < 7 < 4/3: The Collins solution, C(-l/9). 

• 4/3 < 7 < 3/2: The Apostolopoulos solution, 7£+(-l/9). 

• 3/2 < 7 < 2: "Extremely tilted" vacuum plane waves, £^,(—1/9) 

The stability of these points for 7 < 3/2 follows from the eigenvalues of the linearised matrix. For 3/2 < 7 < 2 
several zero-eigenvalues occur and a centre manifold analysis is needed to determine the late-time behaviour. 

4.1.1. The case 3/2 < 7 < 2: the centre manifold. Let us present the centre manifold analysis of the 
equilibrium point £^(—1/9) in some detail. The centre manifold in this case is a 2-dimensional submanifold 
of the 5-dimensional extremely tilted invariant subspace N + (VI—i/g)\ v _ 1 . To find the centre manifold, we 
will therefore set V = 1. Let us choose variables 



(4.1) 



(E + ,E 2 3,^, A, v 2 ) 



3 %/3 V5 

- - + Xt , x 2 , — + x 3 , x 4 , — 



x 5 



let £_, £12, £13 and fl be determined from the constraint equations, and v\ = 1 — v\ ~ v 3- Let us define 
the column vector x = [x\, X2, %3, X4, X5] T . We can now expand the equations of motion to 2nd order in x: 

(4.2) x' = Ax + C(x,x)+0(x 3 ), 

where C(— ,— ) is a bilinear vector-valued function. It is convenient to align the vector x with the Jordan 
canonical form of A. This can be accomplished by defining 



(4.3) 



P = 



7 

5 


2 
5 








■ 





1 








1 




2V3 











Ay/3 




-4V3 





4%/3 


-4\/3 


25 


64 V5 
45 


112V5 

225 









Now, defining y = P : x, the equation for y becomes 

y' = Jy + C(y,y) 

where J is the Jordan block matrix 

J = diag [0, — 



,0, 



The centre manifold correspond to the zero-eigenvalues of J. We can therefore parameterise the centre 
manifold using the variables (2/1,2/4). The next step is to expand the variables yi, i = 2,3,5, in terms of 
(2/1, 2/4) on the centre manifold, to second order. We therefore define the quadratic forms ^(2/1,2/4), i = 2,3, 5, 
such that yi — Yi (2/1,2/4) — is (to second order) an invariant subspace. On the centre manifold we then 
have: 

(4.4) y i = Y i (y 1 ,y±)+0(y 3 ). 

By substituting these into the 2/1 and 2/4 evolution equations, we finally get, on the centre manifold: 

y[ = -162/? + 0(y 3 ), 

(4.5) 2/4 = -162/12/4 + G(y 3 ). 
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To lowest order, these equations can be solved to give y\ w l/(16r), 7/4 w C/r, where a constant of 
integration has been eliminated by a translation of time. Note that, in principle, there are essentially two 
kinds of behaviour: for r < the variables diverge, while for r > the variables decay. However, requiring 
il > leaves only the decaying mode as physically acceptable. These decay rates will therefore be the 
dominant ones since the centre manifold will dominate the behaviour at late times. 

In terms of the original variables, we therefore get the decay rates in N + (VI_i/g) (only dominant decay 
rates are included): 



4 V 60r J 



\/3 / 



1 



1 



4 V 60r 
15 



-42 



43 



J 23 



N 



60r ' 
(1 + 16C)V5 
16072 ' 
(1 + 16C) 
8t 2 ' 

4 V 10r 
(1 + 16C)V3 



4r 



3 

407' 



2 / 9 
vi w — H 

3 V 20r 

V5 / 9 
W2 " -I 1 - 257 



(4.6) Vl-^ 2 « (v/l-V^oexp 



5(2 7 -3) ' 
3(2-7) T 



This confirms that the 'extremely tilted' vacuum plane wave £^(—1/9) is the attractor for N + (VI^i/ 9 ), 
3/2 < 7 < 2. Note, however, that this solution is unstable in the fully tilted space T(VI_i/ 9 ), which can be 
verified by calculating, for example, D' w (3/4)1) close to £^(—1/9). 

4.2. The invariant subspace N~(VI_i/q): non-vortic type VI_i/g models. 

• 2/3 < 7 < 10/9: The Collins solution, C(-l/9). 

• 7 = 10/9: The Wainwright solution, W. 

• 10/9 < 7 < 2: The Collinson-French solution, C^ 7 . We have the following refinement for the 
asymptotic tilt (which follows from an analysis of the eigenvalues): 

★ 10/9 < 7 < (24 - \/6)/15: The tilt is asymptotically zero (CTo). 

★ (24 - v/6)/15 < 7 < 14/9: The tilt is either asymptotically zero (C^o) or extreme {ECTS). 

* 14/9 < 7 < (24 + a/6)/15: The tilt is either asymptotically intermediate {CTi) or extreme 

* (24 + v / 6)/15 < 7 < 2: The tilt is asymptotically extreme (££F_). 

4.3. r(l / /_i/ 9 ): The fully tilted type VI_ 1 / g model. For the fully tilted model, the late-time behaviour 
can be summarised as follows: 

• 2/3 < 7 < 10/9: The Collins solution, C(-l/9). 

• 7 = 10/9: The Wainwright solution, W. 

• 10/9 < 7 < 2: The Collinson-French solution, OF. We have the following refinement for the 
asymptotic tilt: 
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* 10/9 < 7 < (24 - v/6)/15: The tilt is asymptotically zero (CjF ). 

★ (24 — v / 6)/15 < 7 < 4/3: The tilt is either asymptotically zero (CFo) or extreme (fCT 7 -). 

★ 4/3 < 7 < 5/6 + V721/42: The tilt is either asymptotically intermediate (CT\+) or extreme 

■k 5/6 + V721/42 < -f < -f H fa 1.47392: The tilt is either asymptotically oscillatory (closed curve) 
or extreme (EOF '_). 

* 1h < 7 < 2: The tilt is asymptotically extreme {EOF-). 

Most of these results comes from an analysis of the eigenvalues of the various equilibrium points. However, 
for 7 = 5/6 + V721/42 the equilibrium point CT\-y undergoes a Hopf-bifurcation and a stable closed curve 
results. The analysis of this Hopf-bifurcation is given below. Interestingly this closed curve co-exists with 
an extremely tilted attractor. 

4.3.1. The Hopf-bifurcation. Our main aim for this section is to show the following. 

Theorem 4.1. There exists a 70 such that for | + ^p- < 7 < 70 there exists a closed orbit, c(r), acting as 
the future attractor for a set of non-zero measure of tilted Bianchi type VI_i/g models. 

To prove this theorem we will first show the existence of a closed period orbit acting as an attractor in a 
particular subset. We consider the invariant subset given by the Collinson- French solution with 2 tilts. We 
introduce (X, Y) = (vi,v 2 ): 

(4.7) Y' = 2^T+^-X^Y. 

We set X — — ^3 2^ 1 ^0 = 5< ~ 3 ~2{3-yy> 2 ~^ > an< ^ P er f° rm the transformation x = X — Xq, y = Y — Y with 
respect to the equilibrium point CF\ + . 

The normal form of a Hopf-bifurcation can be written 

(4.8) Z 1 = (\ + b\Z\ 2 )Z, 

where b is some complex number and A = a + if} is a parameter. If Re(6) < for a = 0, there exists a stable 
closed orbit for < a sufficently small. 

We will therefore set 7 = | + and expand to cubic terms in x and y. It is also convenient to introduce 
a complex variable z chosen such that it aligns with the Jordan form of the linearised matrix. This can be 
achieved by setting: 

z = x + iay, 

where a is a real number chosen such that the linear term of z' = f(z, z) is 



(4.9) d z f(0 1 0) = ^-1205 + 45^721, 9,-/(0,0) = 0. 
By a transformation 

Z = z + anz 2 + ai 2 zz + a 2 iz 2 + a m z 3 + an 2 z 2 z + a^zz 2 + a 2 22^ 3 , 
we can choose the coefficients o^- and a^-fc so that the equation for Z takes the form 

(4.10) Z' = (\ + b\Z\ 2 )Z + (D(\Z\ 4 ), 
where 



A = - \/ -1205 + 45V72T, Re(b) = (-15990233 + 595193^721) « -0.029. 



Hence, there exists a 70 such that there exists a closed stable orbit for | + -^ji < 7 < 70. 
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The next step is to show that this orbit is also stable in the fully tilted type VI_i/g models. We can show 
this as follows: for a function B, we introduce the average, {B), with respect to the closed orbit c(r), defined 

by 

(B) = i / Bdr, T=l dr. 



T Jc(t) Jc(t) 

We can use this average, using similar manipulations as for the closed curves in the type IV, Vlh and VII/j 
models [5J [TT] to show that: 

Theorem 4.2. Assume that there exists a closed periodic orbit c(r) for the dynamical system ^4.7\ ). Then 

V6(3 7 -4) 5(57-6) 
w 2(3-7) 3(3- 7 ) 

Proof. From the Y equation, we get (T) = — 2^/6 (X) /3. A manipulation of the V equation yields (T) = (S), 
and (S) = (37 — 4) — 2(7 — 1)A (X). These can now be solved to yield the desired value for (X). A similar 
manipulation of the Q equation yields (An). D 

This implies that the vacuum solution is stable when the closed curve is perturbed by fi. Hence, since 
the Collinson-French solution is stable with respect to vacuum perturbations, this closed curve is stable for 
the fully tilted type VT_i/g models whenever | + -^21 < 7 < 7 . 

It remains to determine the maximal value for 70? It seems that this limiting value is related to the 
existence of a heteroclinic orbit originating and ending at the saddle points CTq and &F%, respectively. For 
the limiting value of 7, which we will call 7^, there exists such a heteroclinic orbit, while for values 7 7^ jh 
no such heteroclinic orbit exists connecting these two equilibrium points. We can use this to numerically 
estimate the value for 7^. Our estimate gives: 

1.473920 < lH < 1.473921. 

Since 5/6 + V721/42 = 1.472653409... this means that the region in which the closed orbit exists is extremely 
small0 

Note also that this implies that these closed orbits co-exist with the extremely tilted attractor EOF-, 
which makes this closed orbit even more difficult to detect. In addition, it appears as if the curves asymptote 
to the closed curve relatively slowly, which implies that the numerics have to run for a relatively long time 
in order to see the late-time asymptote. A numerical plot of some generally tilted type VL^j models 
approaching this closed orbit is shown in Fig[TJ 

5. Conclusion 

In this paper we have examined the tilted Bianchi type VI_iyg model in some detail. This model is a 
special case (due to the vanishing of one of the constraint equations) and necessitates a separate analysis 
from the general type VI/j models. We showed that in these models there exists a tiny region of parameter 
space where there exists a closed curve acting as an attractor. This closed curve co-exists with an extremely 
tilted attractor. In the case of the most general irrotational models this closed curve, which appears in terms 
of the tilt velocities, is absent. We have also confirmed the analytical results with an extensive numerical 
investigation. 



However, this seems to be typical for these models, the loophole for the type IV and Vllfo models also appears to be 
extremely small. 
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